Asymmetric polarity reversals, bimodal field distribution, and coherence resonance in 
a spherically symmetric mean-field dynamo model 
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Using a mean-field dynamo model with a spherically symmetric helical turbulence parameter a 
which is dynamically quenched and disturbed by additional noise, the basic features of geomagnetic 
polarity reversals are shown to be generic consequences of the dynamo action in the vicinity of 
exceptional points of the spectrum. This simple paradigmatic model yields long periods of constant 
polarity which are interrupted by self-accelerating field decays leading to asymmetric polarity rever- 
sals. It shows the recently discovered bimodal field distribution, and it gives a natural explanation 
of the correlation between polarity persistence time and field strength. In addition, we find typical 
features of coherence resonance in the dependence of the persistence time on the noise. 

PACS numbers: 47.65.+a, 91.25.-r 



The Earth's magnetic field is known to undergo irreg- 
ular polarity reversals, with a mean reversal rate that 
varies from zero in the Permian and Cretaceous super- 
crons to (4-5) per Myr in the present Q. Typically, these 
reversals have an asymmetric (saw-toothed) shape, i.e. 
the field of one polarity decays slowly and recovers very 
rapidly with the opposite polarity, possibly to rather high 
intensities Q- A general correlation between the persis- 
tence time and the field intensity has also been suspected 
since long [j| . A recent observation concerns the bimodal 
distribution of the Earth dipole moment with two peaks 
at about 4 x 10 22 Am 2 and at about twice that value 0. 

The explanation of these phenomena represents a great 
challenge for dynamo theory and numerics. Remarkably, 
the last decade has seen three-dimensional numerical sim- 
ulations of the geodynamo with sudden polarity reversals 
as one of the most impressive results (cf. for a recent 
overview) . 

Despite the fact that those simulations exhibit many 
features of the Earth's magnetic field quite well, and be- 
cause they do so in parameter regions (in particular for 
the Ekman and the magnetic Prandtl number) that are 
far away from the real ones, there is a complimentary 
tradition to identify the essential ingredients of reversals 
within the framework of simplified dynamo models. This 
has been done, e.g., in the tradition of the celebrated 
Rikitake dynamo model of two coupled disk dynamos Q . 
Another approach has been pursued by Hoyng and col- 
laborators |j| who studied a prototype nonlinear mean- 
field dynamo model which is reduced to an equation sys- 
tem for the amplitudes of the non-periodic axisymmet- 
ric dipole mode and for one periodic overtone under the 
influence of stochastic forcing. Interestingly, even this 
simple model shows sudden reversals and a Poissonian 
distribution of the polarity persistence time. However, 
an essential ingredient of this model to account for the 
correct reversal duration and persistence time is the use 
of a large turbulent resistivity which is hardly justified, 
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FIG. 1: Magnetic field evolution for vanishing noise (D — 0) 
and different values of C. "Field" corresponds to si(r — 0.95). 



at least not by the recent dynamo experiments |8j . 

Sarson and Jones [!j had noticed the importance of the 
transition from non-oscillatory to oscillatory states for re- 
versals to occur. It is our goal to understand this process 
in more detail by studying a simple mean-field dynamo 
model. We focus on the magnetic field dynamics in the 
vicinity of "exceptional points" [l(| of the spectrum of 
a non-selfadjoint operator. We will show that the main 
characteristics of Earth magnetic field reversals can be 
attributed to the square-root character of the spectrum 
in the vicinity of such exceptional points, where two non- 
oscillatory eigenmodes coalesce and continue as an oscil- 
latory eigenmode. 

Our starting point is the well known induction equa- 
tion B = V x (aB) + (/iof) -1 AB for a mean-field dy- 
namo with a helical turbulence parameter a [ll|. acting 
in a fluid with electrical conductivity a within a sphere 
of radius R. The magnetic field has to be divergence- 
free^ ■ B — 0. Henceforth, we will measure the length 
in units of R, the time in units of hqijR 2 , and the param- 
eter a in units of (/lofffl)" 1 . Note that for the Earth we 
get a time scale fiocrR 2 ~ 200 Kyr, giving a free decay 
time of 20 Kyr for the dipole field. 

We decompose B into a poloidal and a toroidal part, 
B = — V x (r x \7S) — r x VT. The defining scalars 
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FIG. 2: Magnetic field amplitude and frequency for D = in 
dependence on C. Note the phase transition at Ct = 1.27893. 



S and T are expanded in spherical harmonics of degree I 
and order m with the expansion coefficients s; m (r, t) and 
tl,m(r,t). In order to allow for very long simulations (to 
get a good statistics for the persistence time and the field 
amplitudes), we consider an a 2 dynamo with a radially 
symmetric helical turbulence parameter a. In [12J we had 
shown that such simple a 2 dynamos can exhibit oscilla- 
tory behaviour in case that a(r) changes its sign a long 
the radius, which is not unrealistic for the Earth |13| . 
For spherically symmetric and isotropic a, the induction 
equation decouples for each I and m into the following 
pair of equations: 
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These equations are independent of the order m, hence 
we have skipped it in the index of s and t. The boundary 
conditions are dsi/dr\ r= \ + (I + l)sj(l) = i;(l) = 0. 

In the following we restrict ourselves to the dipole field 
with I = 1 (the influence of higher multipole fields will be 
considered elsewhere). At first, we choose a particular ra- 
dial profile a(r) which can yield an oscillatory behaviour 
[l2^ . Magnetic field saturation is ensured by quenching 
the parameter a(r, t) with the angular averaged dipole 
field energy which can be expressed in terms of sx(r, t) 
and ti(r, t). In addition to that, we perturb the a-profile 
by noisy "blobs" which are assumed constant within a 
correlation time r. In summary, a(r, t) takes the form: 



a(r,t) = C- 



-21.5 + 426.4 r 2 - 806.7 r J + 392.3 r 4 



l + E 



d(r Sl {r,t)) 
dr 



t 2 (r,t) 



(3) 



where the noise correlation is given by + = 

D 2 {\ - |ii|/r)e(l - |fi|/T)<Jy. C is a normalized dy- 
namo number measuring the overcriticality, D is the 
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FIG. 3: (a) Details of the field evolution for C = 1.2789 (cf. 
Fig. 1), with growth rates and frequencies (inset) resulting 
from the instantaneous, quenched a profiles, (b) Growth rate 
curves for instantaneous a profiles scaled by C' . "K" denotes 
the kinematic (unquenched) a profile with the (encircled) crit- 
ical point at C' — 1. 



noise amplitude, and E is a constant measuring the in- 
verse equipartition field energy. 

The equation system (l)-(3) is time-stepped using an 
Adams-Bashforth method with a radial grid spacing of 
0.02 and a time step length of 2 x 10~ 5 . For the following 
examples, the correlation time r has been set to 0.02, and 
E has been chosen to be 0.01. 

For the noise- free case (D = 0), we show in Fig. 1 
the evolution of the magnetic field for different values of 
C (hereafter, "field" stands always for the value of s\ 
at r = 0.95). At the slightly overcritical value C = 1.1 
we observe a nearly harmonic oscillation with a frequency 
close to one. With increasing C, this frequency decreases, 
and at the same time the signal becomes saw-toothed (or 
better: "shark-fin shaped") and eventually rectangular. 
At the critical point Ct = 1.27893 a phase transition 
to a non-oscillatory dynamo occurs, which can also be 
identified in the frequencies and field amplitudes shown 
in Fig. 2. 
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FIG. 4: Illustration of the various phases of reversals in the 
vicinity of an exceptional point of the spectrum. In the 
noise free case (D=0) oscillatory and steady dynamos can be 
sharply distinguished. For D > 0, we get a bimodal behaviour 
around the maximum of the non-oscillatory branch. 



In Fig. 3, we analyze the field evolution for C = 1.2789 
in detail. Figure 3a shows the field during one reversal, 
together with the growth rate and the frequency that re- 
sult from the instantaneous, quenched a-profile. Eight 
of these profiles (at the moments 1-8), together with the 
unquenchcd profile, are scaled by a factor C yielding 
the growth rate curves of Fig. 3b. These curves may 
help us to identify the following phases for a reversal: 
a slowly starting, but self-accelerating field decay in the 
non-oscillatory branch (points 1 and 2), followed by the 
actual polarity reversal in the oscillatory branch (between 
points 3 and 4), a fast increase of the field (points 4-6) 
which results in a rapid a-quenching. This rapid quench- 
ing drives the system back to the point 8 which basically 
corresponds to the initial point 1, only with the opposite 
field polarity. A peculiarity of our particular model is the 
existence of a second exceptional point beyond which the 
dynamo becomes non-oscillatory again (between points 
5 and 6), which is, however, not crucial for the indicated 
reversal process. 

In an attempt to simplify this reversal picture, we con- 
sider in Fig. 4 all growth rate curves of Fig. 3b as being 
collapsed into a single one (for that purpose we neglect 
their slight shape differences and consider only a shift). 
By this "optical trick", the various phases of reversals 
can be visualized as a "move" of the actual growth rate 
point relative to this collapsed growth rate curve. Figure 
4 explains also the phase transition between oscillatory 
and steady dynamos and the discontinuity of the field 
amplitude from Fig. 2. The local maximum of the non- 
oscillatory branch of the growth rate curve is unstable 
and repels the dynamo in two different directions, de- 
pending on whether the growth rate there is negative or 
positive. If it is negative, the field decays, we get weaker 
a quenching, and a subsequent reversal. If it is positive, 
the field increases, we get stronger a quenching, hence 



FIG. 5: Time series for C = 1.27 and C = 1.29 at D = 0.5. 
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FIG. 6: Field histogram for C = 1.27 and C = 1.29. 



the system moves to the left of the maximum. 

What happens now if we switch on the noise? The 
influence of the noise is quite different for values of C 
below and above Ct- For C > Ct, the noise is the 
only possibility to trigger a move from left of the local 
maximum to the right, hence the persistence time will 
decrease from infinity to some finite value. For C < Ct, 
the noise will sometimes push the maximum above the 
zero growth rate line allowing the system to jump to the 
left of the maximum where it can stay for a while. Hence, 
we will get an increase of the persistence time. 

For a moderate noise intensity D = 0.5, we depict in 
Fig. 5 the time evolution for the values C = 1.27 and 
C = 1.29 which are slightly below and above Ct, respec- 
tively. First of all, a drastic difference between the per- 
sistence times is still visible. The duration of a reversal, 
however, is quite identical for both values of C, although 
its sensible definition is not obvious. With a side view on 
the real geomagnetic field (for which the dipole compo- 
nent should decay approximately to the strength of the 
non-dipolc components before a reversal can be identi- 
fied), we tentatively define the reversal duration as the 
time the field spends between ±1/4 or ±1/10 of its mean 
intensity. In the first case we get a reversal duration of 
30 Kyr, in the second case of 12 Kyr, which comes close 
to the observational facts. 
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FIG. 7: Persistence time (T) and its normalized standard 
deviation a(T)/{T) for different C in dependence on D. For 
each point, a minimum of 10 4 reversals was used. 



get typical persistence times of the order of 200 Kyr (and 
larger), and a typical reversal duration of 10-30 Kyr. 

We point out that this model is not a physical model 
of the Earth dynamo, in particular owing to the miss- 
ing North-South asymmetry of a. Nevertheless, we are 
convinced that the generic features of dynamos in the 
vicinity of exceptional points as they were illustrated in 
this paper can also be identified in much more elaborated 
dynamo models. 
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For both values, C = 1.27 and C = 1.29, a bimodal 
behaviour of the dynamo can be observed in the field 
histograms of Fig. 6. For C = 1.27 and D = 0.15, 
we observe a clear double peak on both polarity sides, 
centered approximately at ±0.18 and ±0.27. Evidently, 
the dynamo is mostly in the weak field state (right of 
the maximum, cp. Fig. 4), with some "excursions" to 
the strong field state (left of the maximum). With in- 
creasing D this double peak is smeared out, leaving a 
very broad maximum for D = 0.5. The three histograms 
for C — 1.29 show a maximum at the strong field value 
and (although not a maximum) a pronounced flatten- 
ing at the weak field value. This means, the dynamo is 
mostly in the strong field state, with some intermediate 
stopovers in the weak field state, from where it can start 
a reversal. 

An overview about the mean persistence time (T) and 
its normalized standard deviation o~(T)/(T) in depen- 
dence on C and D is given in Fig. 7. Note the drastic 
difference that small values of noise have on the dynamo 
behaviour for C < Ct and C > Ct- The normalized 
standard deviation, at least for the curves with C > 1.26, 
has a clear minimum around D — 3 which represents a 
typical feature of coherence resonance |l4j . We only note 
here that recently stochastic resonance models have been 
discussed with view on a possible triggering of field rever- 
sals by the 100 Kyr period of the Earth orbit eccentricity 

To summarize, our simple spherically symmetric a 2 
dynamo model shows that asymmetric polarity reversals 
and bimodal field distributions are generic features of the 
dynamo behaviour in the vicinity of exceptional points. 
Using only the molecular resistivity of the outer core, we 
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